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INTRODUCTION 


Tli  e orient  at  ion  of  t h o s cns  i n p axis  o I a n I its  t r u in  out  wit  h 
res  poet  to  a 1 ixed  coordinate  system  in  space  is  ret  erred  to 
as  its  attitude.  From  the  attitude  of  a sensor,  one  may 
determine  the  angle  it  makes  with  any  other  vector  such  as 
the  rocket's  velocity  vector,  the  sun  line  vector,  etc.  For 
certain  sensors,  meaningful  interpretation  of  the  instrument's 
output  can  be  accomplished  only  when  the  attitude  of  the  sensor 
ax  is  is  known . 

The  first  section  of  this  report  discusses  various  types 
°f  gyroscopic  platforms  used  as  attitude  measuring  systems  on 
vehicles.  The  systems  include  the  Attitude  Control  System  with 
two  free  gyros,  the  Miniature  Attitude  Reference  System  with 
segmented  lengths  over  a 5 volt  span  for  the  roll,  pitch  and 
yaw  axes,  and  the  Miniature  Inertial  Digital  Attitude  System 
with  a one-to-one  correspondence  between  angular  displacements 
and  digital  coding  outputs. 

The  next  discussion  centers  upon  processing  procedures 
which  produce  final  attitude  information  from  vehicle  gyro 
measurements.^  An  analvsis  is  shown  which  relates  any  onboard 
probe  vector  P to  the  local  coordinate  system  fixed  at  the 
launch  site.  With  the  orientation  of  P determined,  the  attack 
angle  with  any  other  attitude  determined  or  predetermined 
vector  can  be  supplied. 

To  provide  continuous  final  output  for  the  attitude  system 
measurements,  data  refinements  techniques  are  introduced  with 
statistical  error  bound  values  on  the  predicted  versus  measured 
output.  Tolerance  levels  are  set  on  the  measured  input  data 
to  ensure  that  all  predicted  data  will  be  within  a specified  a 
deviation  . 

Oftentimes  the  initial  orientations  of  onboard  sensors 
require  corroboration.  A phase  and  time  shift  study  is  intro- 
duced in  the  third  section  to  compare  the  phase  relationship 
between  measured  and  attitude  predicted  magnetometer  pitch 
angles  during  well-behaved  areas  of  a vehicle's  flight.  The 
routines  can  be  used  for  both  side  and  axial  magnetometer 
output . 

Finally,  some  of  the  more  common  types  of  problems  with 
gyro  systems  are  discussed.  The  correction  procedures  used 
on  these  malfunctions  are  discussed  and  the  results  exhibited. 
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SUCTION  1 


GYROSCOPIC  PLATFORMS  USLD  AS 
ATTITUDE  MEASURING  SYSTEMS 

1.1  Attitude  Control  System 

Space  Vector  Corporation's  Attitude  Control  System  (ACS)  is  comprised  of 
a Roll  Stabilized  Platform  (RSP)  with  two  free  gyros  which  output  a yaw,  pitch 
and  roll  vehicle  motion  from  top  dead  center  (TDC)  of  the  gyro.  TDC  is  the 
reference  for  all  onboard  probes  (Fig.  1).  The  roll  and  yaw  are  true  measure- 
ments of  vehicle  motion  but  since  there  are  only  two  free  gyros,  true  vehicle 
pitch  is  a function  of  yaw.  Further,  if  the  vehicle  yaws  over  85°,  this  will 
cause  loss  of  initial  reference  for  the  system  and  make  attitude  determination 
almost  impossible. 


TDC  Reference 


Figure  1 

Various  models  of  the  ACS  were  used  as  attitude  measuring  systems  on 
vehicles  analyzed.  They  all  function  similarly  but  output  differently.  Typical 
examples  are  given  below. 
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Model  10580  ACS  flown  on  IC(>50.02-1A  output  pitch  course,  pitch  fine,  yuw 
coarse,  yaw  fine,  and  roll  gyro  data.  The  pitch  and  yaw  coarse  data  were 
calibrated  in  .12  volt  increments  with  a linear  conversion  range  of  190° 
corresponding  to  a 5 volt  span.  The  pitch  and  yaw  fine  data  over  a similar 
volt  span  had  a conversion  range  of  ±5°.  The  roll  telemetry  data  had  a linear 
conversion  range  of  ±180°  on  a 4.8  volt  span. 

According  to  Figure  2,  pitch  course  for  this  model  was  decreasing  voltage 
output  when  pitching  down,  pitch  fine  was  increasing  voltage  when  pitching 
down.  A yaw  right  motion  was  increasing  yaw  coarse  voltage  and  a decreasing 
yaw  fine  voltage.  A clockwise  roll  viewed  aft  to  nose  was  increasing  voltage. 

Model  10390  ACS  flown  on  1C519.07-1B  output  only  pitch  coarse,  yaw  coarse 
and  roll  gyro  data.  The  pitch  data  was  calibrated  over  a -+90°  corresponding 
to  a 5 volt  span.  However,  to  model  the  calibration  information,  discrete 
polynomial  expressions  were  used  over  finite  intervals.  From  0 volts  - 4 volts 
a linear  model  was  used  and  from  4 volts  - 5 volts,  a second  degree  model  best 
represented  the  calibration  information.  The  yaw  data  was  converted  to  angle 
measurements  in  a fashion  similar  to  that  used  for  the  coarse  yaw  output  on 
model  10380.  The  roll  telemetry  data  had  a linear  conversion  range  of  four 
t45°  segments,  each  on  a 5 volt  span. 

Pitch  and  yaw  coarse  output  reterences  were  s'imilar  to  those  for  model 
10580,  but  clockwise  roll  output  was  decreasing  voltage  when  viewed  aft  to  nose 

Modei  10550  ACS  flown  on  A10. 304-1  output  pitch  coarse,  pitch  fine,  yaw 
coarse,  yaw  fine  and  roll  gyro  data.  The  pitch  coarse  segment  was  calibrated 
from  -20°  to  +170°  over  a 5 volt  span  with  decreasing  voltage  when  pitching 
down.  The  yaw  coarse  segment  was  calibrated  for  a range  of  ±70°  over  a 5.3 
volt  span,  i.e.,  .8  to  4.1  volts.  A yaw  right  motion  was  decreasing  voltage 
for  the  coarse  yaw  output.  Three  data  points  were  recorded  for  pitch  fine 
data,  and  the  curve  connecting  each  pair  of  points  is  essentially  a straight 
line.  These  points  are: 
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ACS  COORDINATE  REFERENCE 


ROLL  ANGLE  INCREASING  CW^ 


YAW  RIGHT 


‘LAUNCH  RAIL  * Odegrees  ROLL 

Figure  2 
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-1  0.41 

0 3.  JO 

+ 1 4.-32 

No  data  calibration  points  were  recorded  for  yaw  fine.  According  to  Figure  2, 
pitch  tine  was  increasing  voltage  when  pitching  down  and  yaw  fine  was  increasing 
voltage  when  yaw  to  the  right.  The  roll  segment  was  calibrated  from  180°  CCW  to 
1,6  Ci\  >ver  a 4./  volt  span  with  voltage  decreasing  for  a CW  roll  wiien  viewed 
aft  to  nose. 


Another  ACS  gyro  reference  system  was  the  Space  Data  Corporation  system  used 
tor  the  vehicle  EX531.43-1.  The  pitch  and  yaw  channel  output  were  ±180°  over  a 
5 volt  span  and  the  roll  channel  output  was  a 360°  segment  over  a 5 volt  span. 
Although  not  used  in  the  data  reduction,  tnc  pitch  and  yaw  high  resolution  were  +10° 
The  zero  reference  for  the  gyro  roll  was  the  launch  rail  (Figure  3).  Gyro  pitch 
zero  was  73.3°  launcher  elevation  setting  at  gyro  uncage  and  gyro  yaw  zero  was 
46.2°  launcher  azimuth  setting  at  gyro  uncage.  Looking  forward,  yaw  right  is 
positive,  pitch  up  is  negative  and  roll  clockwise  is  positive. 


ka  1 1 

o 


180 

Figure  3 


1.2  MARS  Attitude  System 

I he  Space  Vector  Corporation's  Miniature  Attitude  Reference  System  (MARS)  pro- 
vided gyro.'copic  output  for  the  determination  of  rocket  attitude  for  the  vehicle 
A10. 302-1.  Hie  roll  axis  coincided  with  the  rocket  axis  at  launch  and  the  orienta- 

i 

tion  of  a gyro  reference  notch  was  used  to  determine  the  location  of  the  yaw  axis. 
The  cross  product  of  the  yaw  axis  with  the  roll  axis  yielded  the  pitch  axis.  The 
axes  and  rotations  of  the  MARS  system  are  identified  in  Figure  4. 
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Top 
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Figure  4 


The  segmented  lengths  over  a 5 volt  span  for  each  MARS  axis  are: 

a)  Roll:  90° 

b)  Yaw:  60° 

c)  Pitch:  45°. 

Segment  identification  of  the  pitch  and  yaw  outputs  for  this  system  requires 
continuous  monitoring  of  the  output  to  determine  the  instance  of  instantaneous 
transitions  from  one  segment  to  the  next.  By  keeping  track  of  the  transitions 
and  their  direction  evidenced  by  a plus  or  minus  5 volt  step  change,  one  may 
correctly  identify  the  segments  being  presented.  The  roll  output  identification 
was  a 5°  *3°  shorted  segment  in  the  center  of  the  1st  segment  which  appeared 
as  a constant  output  voltage  for  that  portion  of  the  segment.  This  short  was 
repeated  every  4th  segment. 


1.3  MIDAS  Attitude  System 

Space  Vector  Corporation's  Miniature  Inertial  Digital  Attitude  System  (MIDAS) 
provided  gyroscopic  output  used  to  calculate  the  attitude  of  many  rockets  such 
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as  AID. 403 


tl 


system  but  the 
S volt  span  rep 
MIDAS  system's 
and  result  in  a 
digital  codings 
in  degrees  is 


lhe  roll,  pitch  and  yaw  references  are  the  same  as  for  a MAKS 
outputs  differ  considerably . Whereas  a MAKS  system  outputs  in  a 
resenting  some  degree  segment  having  a predetermined  length,  the 
roll,  pitch  and  yaw  measurements  are  digitized  by  optical  encoders 
one-to-one  correspondence  between  angular  displacements  and 
. for  the  digital  coding  variable  n,  the  angular  displacement  <? 


1024 


x 360 


No  bias  values  are  needed  for  conversion  since  the  MIDAS  system  represents  the 
displacement  from  the  uncaged  position  of  the  gyro  which  fixes  the  coordinate 
reference  system.  The  flow  of  attitude  data  for  MIDAS  vetucles  is  simplified 
considerably  when  compared  with  the  flow  of  ACS  and  MARS  vehicles. 


SECTION  2 


PROCESSING  PROCEDURES,  MATHEMATICAL  ANALYSIS,  AND 
REFINEMENT  TEGHNI-QUES 


2.1  Processing  Procedures  for  all  Rockets 

The  oscillogram  containing  vehicle  gyro  measurements  is  inspected  for 
calibrations,  extreme  noise,  gyro  malfunctions,  approximate  spin  rates,  pro- 
cessional period  and  half  cone  angle.  The  digital  tape  containing  the  gyro 
data  is  then  unpacked  and  converted  to  degrees  and  plotted.  These  plots  are 
then  compared  with  results  of  the  oscillogram  inspection. 

If  no  gyro  malfunction  occurs  then  the  quick-look  attitude  and  angles  of 
attack  for  the  rocket  axis  are  calculated  based  upon  preliminary  information. 

If  the  attitude  measuring  system  is  a MIDAS,  then  attack  angles  for  requested 
side  probes  are  calculated  at  the  same  time  as  attack  angles  for  the  rocket 
axis.  If  attitude  for  side  probes  is  requested  for  an  ACS  or  MARS  systt..., 
the  roll  data  is  converted  to  degrees,  checked  for  proper  conversion  and  then  th 
preliminary  side  probe  attitude  is  calculated. 

All  extreme  noise  and  calibrations  are  filtered  from  the  gyro  conversions 
and  this  output,  if  displaying  a lack  of  continuity,  is  curve  fitted  by  tech- 
niques discussed  in  Section  2.3. 

Finally,  if  onboard  probes  measure  magnetic  field,  lunar,  or  solar  intensi- 
ties, it  is  possible  to  corroborate  the  calculated  attitude.  Depending  upon 
the  accuracy  of  the  onboard  measurements,  an  error  bar  estimate  of  the  attitude 
can  be  included  as  part  of  the  normal  output.  A flow  of  the  attitude  data  for 
both  ACS  and  MARS  vehicles  is  displayed  in  Figure  5,  and  a flow  of  the  data 
for  MIDAS  vehicles  in  Figure  b. 

2.2  Analysis 

^ A A 

Given  the  unit  vectors  X,  Y,  Z in  the  directions  of  the  gyro  roll, 
and  yaw  axes,  the  direction  of  the  longitudinal  axis  of  the  vehicle  can 
pressed  as  the  vector  e": 


pitch 
be  ex- 
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L-  LOW  01-  ATT  ITU Dl:  DATA  1 OR  MIDAS  VEHICLES 


4 

(1)  = Raw  Data  Tape 

(2)  “ Oscilloyram 

(3)  = Information  Sheet 

(E)  = Trajectory  Tape 


Figure  6 
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Z cosy  snip 


(2.1J 


e"  = X cosy  cops  + Y siny  +■ 


where  y and  p are  the  true  vehicle  vaw  and  pitch  respectively.  Now  let  e,,  , 

e-  , er  unit  vectors  in  the  directions  of  true  North,  hast  and  the  local 
vc  c 

vertical,  fixed  at  the  launch  site.  Within  this  system,  the  elevation  of  a 
vector  .nay  be  determined  as  the  angle  it  makes  with  the  local  horizontal  plane 
with  the  direction  of  the  vertical  taken  as  positive  elevation.  The  azimuth 
of  a vector  is  the  angle  between  its  projection  in  the  horizontal  plane  and 
eQc  measured  positive  clockwise  (see  Fig.  7). 


(hast) 


8 = Llevation  of  e" 
r 

<{>  = Azimuth  of  e" 
r 


F igure  7 


Using  the  direction  coefficients  of  the  gyro  axes  at  launch,  the  unit  vectors 

X,Y,Z  can  be  expressed  as  linear  combinations  of  the  eartn-based  system  °tc>ej>c> 

er  . If  we  call  the  coefficient  matrix  11,  then 
* c ’ 


Nc 
= B iC0 


According  to  equation  (6)  of  Reference  [3],  we  can  define  a system  of  orthonormal 
vectors  c” , e j , e' ' at  any  time  in  flight  as  a linear  combination  of  the  gyro  axes. 
If  A represents  the  coefficient  matrix  of  the  gyro  output,  then 
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— - 

c" 

X 

r 

o" 

= A 

Y 

l 

e" 

Z 

2 

In  particular,  as  described  in  (2.1J 


a = cosy  cosp 


a = smy 
12 


a = cosy  sinp  . 

1 3 


To  transform  the  time  dependent  e"  e",e"  system  into  the  local  eB  ,c„ 

r i 2 vc  rc 

system,  form  the  matrix  product  of  A and  B.  Letting  the  matrix  C = AB  then 


c- . = l a.  b . (l, j = 1 ,2,3) 
lj  . lm  mi  J > > j 

J m=l  J 


and 


e" 

r 

e0c 

e" 

l 

u 

II 

e4>c 

e" 

2 

yc_ 

(2. 2) 


Once  the  attitude  ot  the  vehicle's  main  axis  has  been  determined  as 
direction  cosines,  the  elevation  and  azimuth  may  be  derived.  From  (2.2J;  we 
have: 


er  = e0  cos®  cos<*>  + e<ti  cos9  sin!*)  + e sin0 
c 9c  rc 


Therefore: 


0 = sin'1  (c13)  , 


12 


a iid 


„ -1  ,^12 

• tan  (- — ) • 

L11 


rile  eievition  and  azimuth  of  e'.'  and  e"  are  similarly  derived. 

In  order  to  extend  this  approach  to  a vector  I’  lying  in  the  sensing  direc- 

A 

tion  of  a probe  having  any  orientation  on  the  rocket,  P is  represented  as 

P = e”  cosk  + e"  cosu  sink  + c"  sinq  sink  , 

where  A is  the  angle  between  e’^  and  I’,  and  u is  the  angle  between  e"  and  the 

projection  of  P in  the  plane  of  eV  and  e".  furthermore,  if  0 and  <i>  are  the 

1 „ 2 P P 

respective  elevation  and  azimuth  of  P,  it  follows  that 


P = en  cos0  cos$  + e.  cos0  s inf  + c sin0 
ec  P P <J>C  P P rc  P 


(2.3; 


he  can  also  define  P 


,P2, 

p,  - 
0 

■ the 

vector 

Ip,  j 

d,  . 

d , 

l 1 

1 

11 

12 

k 

u22 

p»i 

d51 

d32 

_ 

1 

components  of  P - by 


d13 

V 

c 

d23 

V 

d„ 

e 

r 

c 

the  transformation: 


(2.4) 


where  each  d. . (i,j=l,2,3)  is  a product  of  c.j  and  the  corresponding  direction 
cosine  of  P.  Using  the  same  technique  as  for  c^,  the  coefficients  of  eg^.e^,, 
er^  in  (2.3)  and  (2.4)  are  equated,  readily  yielding  0^  and 


finally,  if  V is  any  vector  associated  with  the  rocket  flight,  the  angle 

/A 

between  P and  V may  be  derived.  The  unit  vector  in  the  direction  tf  V is  given 
by 


c0  V1  + c<p  v2  + er  v3 
c _c c 

"Tvi " 


C0 


+ e. 


+ e 


vi 
r 3 
c 


(2.5) 
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where  the  \r  (.1  = 1, 2,5)  terras  are  the  components  of  V in  the  directions  of  north, 
east  and  vertical;  the  resulting  vl  terras  are  thus  the  direction  cosines  of  V. 
The  scalar  product  of  (2.3)  and  (2.5)  yields 

A A 

l’*V  = v!  cos n cos!)  + v cosO  s i n2>  + vl  sinO 

1 p P 2 p yp  3 p 


Hence,  the  angle  i p,  between  P and  V is  given  by: 


-1 


ii»  = cos  (vl  cos0  cos<l>  + vl  cosO  sin<&  + vl  sinO  ) 
1 P p 2 p vp  j p' 


2.3  Data  Refinement 


Due  to  noise  or  in-flight  calibrations  often  present  in  the  recorded  data, 
filtering  procedures  are  modified  or  developed  as  required  to  provide  continuous 
final  output  for  the  attitude  determination  system.  The  primary  routines  used 
are  Fourier  Series  and  n^h  degree  polynomial  approximations. 


To  provide  smooth,  continuous  and  acceptable  attitude  information  for 
vehicles  when  they  have  established  angular  precessional  velocity  and  a half 
cone  angle  it  is  possible  to  generate  gyro  yaw  and  pitch  information  by  means 
of  the  Fourier  Series  expansion.  During  this  well-behaved  areas  of  a particular 
vehicular  flight,  the  coarse  and/or  fine  yaw  and  pitch  data  can  be  predicted  by 
tne  Fourier  Series 


y 


VJ  V* 

— + z 

n=l 


(a  cos 
n 


117TX 


, Xttx. 

b sin  — ; — ) 
n L 


Tnis  problem  has  proven  successful  due  to  the  periodic  nature  of  the  data. 
However,  for  most  rocket  flights,  quick  convergence  of  the  Fourier  Series  is 
obtained  with  the  approximation 


y = fj  + f^t  + fj  sincot  + f^  coscot 


The  addition  of  the  linear  term  is  needed  to  account  for  an  occasional  linear 
shift  in  the  converted  data.  A preliminary  angular  velocity  co  is  selected  from 
a study  of  the  oscillograms  of  the  raw  data,  and  is  further  refined  by  an  option 
within  the  fitting  routine. 


1 A 


During  regions  in  which  .1  vehicle  is  not  coning,  yaw  and  pitch  outputs  are 
separated  into  discrete  time  intervals.  These  intervals  are  then  curve  fitted 
with  polynomials  up  to  the  20tn  degree  when  necessary.  The  n^1  degree  polynomial 
routine  calculates  an  RMS  value  between  measured  and  predicted  data  for  each  of 
tne  polynomial  approximations.  The  minimum  RMS  value  determines  the  degree  of 
the  polynomial  approximation  to  be  used  in  the  specified  interval. 


The  PMS  is  normally  defined  as 


RMS 


(x.-x)2 

X 


(2.6) 


where  x(  is  the  difference  between  a measured  and  predicted  data  value.  To 
simplify  computation  and  programming,  an  equivalent  form  of  (2.6) 


RMS  = 


_ 1 
x“ 


(2.7) 


was  used,  lor  a reasonably  large  sample  size  N,  the  RMS  deviation  approximates 
the  standard  deviation  a.  Since  a £30  deviation  accounts  for  * 9y.79o  of  tne 
measured  data  in  a normal  distribution,  tolerance  levels  are  set  on  input  data 
to  ensure  that  all  predicted  data  will  be  within  the  specified  0 deviation  while 
utilizing  the  results  of  (2.7). 

figures  (.8-22)  show  normal  angles  of  attack  using  data  refinement  techniques. 
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Slid  I ON  3 


PHASE  A AID  TiML  SHIFT  STUDIES 


1o  corroborate  the  initial  orientations  of  on-board  sensors,  measured 
and  attitude  predicted  magnetometer  pitch  angles  are  compared.  A software 
package  was  developed  to  compare  the  phase  relationship  during  well-behaved 
areas  of  a vehicle's  flight. 

This  package  involves  a least  squares  study  on  the  magnetometer  data 

(refer  to  Section  3.2  for  a detailed  description) . The  output  produced  includes 

a table  of  maximum  (minimum)  times  T measured  data  values  in  degrees, 

associated  AT  values  and  spin  rates  for  the  corrected  magnetometer  outputs. 

Also  included  are  a table  of  maximum  (minimum)  times  T . predicted  data 

gen  1 

values  in  degrees,  associated  AT  values  and  spin  rates  for  the  predicted 

gen 

magnetometer  output  from  the  attitude  data.  The  two-time  dependent  tables  are 
merged  for  all  consecutive  time  values  of  each  array  and  the  phase  shift 
between  the  associated  values  are  printed. 

Tor  T < T + AT  . the  phase  shift  x between  the  raw  and  the  pre- 
gen  — meas  meas’  1 A ^ 

dieted  magnetometer  data  is 


X 


T - T 

gen  meas  0 

— x obO 

AT 

meas 


For  0 < x < lbO°  the  predicted  data  is  lagging  or  coincident  with  the  measured 
magnetometer  data.  For  180°  < X < 360°,  the  measured  data  is  said  to  lag  the  pre- 
dicted by  360°  - x-  The  following  figure  depicts  the  relationships  just  discussed 


T T 

meas  gen 
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In  the  previous  illustration,  the  predicted  data  from  the  attitude  is  lagging 
the  raw  data  since  we  are  assuming  that  (T„en  - Tmcas)  1 ^-ATm  s* 

To  check  for  possible  timing  problems,  a time  shi ft  between  measured  and 
predicted  data  is  computed.  Statistical  studies  include  a mean  time  shift  and 
associated  RMS  value  for  each  specified  time  interval. 

3.1  Bias  liquation 

In  order  to  accurately  determine  the  times  of  maximum  magnetometer  output, 
smooth  magnetometer  readings  are  required.  Normal  data  transmission  does 
not  provide  such  smooth  output,  so  to  allow  for  proper  usage  of  the  routines 
to  be  described  in  Section  3.2,  a curve  fit  was  put  through  the  bias  values 
of  the  magnetometer  data.  The  samples  selected  are  taken  from  the  coning 
portion  of  flight  - i.e.,  from  approximately  40  kins . on  ascent  to  approximately 
40  kms.  on  descent. 

Three  equally  spaced  time  intervals  are  selected  to  obtain  a separate  bias 
approximation  B for  each  interval  from  the  equation 

IV. 

_ l 

n 

A typical  study  will  produce  a value  of  n in  the  order  of  2 * 10  . 

The  standard  bias  equation  as  a function  of  time  is  normally  approximated 
by  a first  degree  polynomial.  However  to  obtain  greater  aqcuraey  with  higher 
order  terms  in  the  fitting  routine,  the  value  of  the  bias  B(tJ  as  a function  of 

> 

B(t j = at*"  + bt  + c . 

IV. 

However,  from  the  discussion  above  we  can  approximate  B(tJ  by  — — in  each  of 
the  intervals  studied.  This,  in  turn,  leads  to  the  system  of  equations 


1 7 


< 

+ btl 

+ c - 

1 

at2 

+ bt2 

+ c = B2 

at  3 

+ bt3 

+ c = B, 

where 


j = i,2,3  is  taken  as  the  mean  time  in  the  jth  interval  and 


B. 

J 


n 

Z 

i = l 


V. 

1 


j=l,2,3 


As  tie  value  of  n increases,  the  maximum  error  that  occurs  in  the  approxi- 
mation of  B^  decreases.  This  fact  can  bo  verified  by  calculations  on  the 
TV. 

on  the  expression  . 

n 


This  procedure  to  determine  the  bias  equation  as  a function  of  time,  forms 
part  of  the  procedure  used  on  all  vehicle  flights  when  examining  phase  shifts 
between  measured  and  predicted  data. 


3.2  Maxima  of  the  Magnetometer  Measurements 

Giver  the  sinusoidal  nature  of  the  lateral  magnetometer  outputs,  the 
prepared  software  routines  have  the  capability  of  determining  either  the  times 
of  relative  maximum  outputs  or  the  times  of  relative  minimum  outputs.  If  a 
least  squares  approximation  of  the  data  is  desired  in  the  maximum  areas  of  the 
output,  then  a test  is  performed  on  the  data  V.  satisfying  the  inequality. 

B(T. ) < V.  . 

1'  1 

Likewise,  if  a fit  of  the  minimum  areas  is  desired,  then  only  those  V values 
satisfying  the  inequality 
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arc  stored  for  computation  ui  the  least  squares  polynomial  fits  of  the  specified 
areas  of  the  data. 

in. order  to  initiate  the  procedure  for  storing  data  in  the  maximum  areas 
of  tiic  least  .uuares  aporoximat  ion.  a value  V must  lie  found  such  that  V it  f t i 

i il 

and  V < b (i  , ) . Once  initiated,  data  is  stored  until  V • lilt  / and 
i i l-l  il 

V(  j B(t ^ where  j > i.  This  last  condition  is  a necessary  one  for  ending 

the  sample  region  hut  not  a sufficient  condition.  The  sufficiency  test  is  made 

bv  examining  V , to  ensure  that  V is  not  a noise  point.  If  V.  , < Bit  .) 

j +1  j * j+1  j+1 

then  a least  sijuarss  parabolic  fit  is  performed  on  the  sampled  data.  However, 
if  B(t.+jJ,  the  point  V ■ is  discarded  and  the  summing  procedure  of  tnc 

least  squares  continues.  Similar  tests  are  made  in  the  minimum  regions  of 
fitted  data  except  that  the  inequality  signs  are  reversed. 

In  each  specified  region,  the  data  is  fit  to  the  quadratic 

Y = At“  + fit  + C . 

Ihe  value  of  the  bias  B(t.J  at  which  the  sampling  begins  and  that  value  B(tij, 

j > i at  which  the  sampling  ends  arc  obviously  just  the  points  of  inflection 

for  the  sinusoidal  curve.  The  time  of  relative  maximum  (minimum)  is  obtained 
dY 

by  setting  = 0 and  solving  for  time  t.  This,  in  turn,  provides  us  with  a 
predicted  output  value 

Y = A'f“  + BTm  + C 


where 


Should  the  measured  and  predicted  outputs  deviate  by  more  than  some 

specified  limit,  usually  taken  to  be  4T,  the  acceptance  lc'vel  for  least  squares 

data  may  tie  modified  so  tiiat  the  sampled  regions  more  closely  approximate  the 

parabolic  fit.  This  modification  amounts  to  accepting  the  area  containing 

maximum  data  when  V • B(t.j  a k where  k is  a constant  determined  'rom  examina- 
i i 

tion  of  the  use 1 1 iograpiis . for  minimum  output  regions,  the  inequality  takes 
the  form  V < B(t.)  - k. 

i v l 


# 


For  each  value  of  generated  by  the  1st  derivative  test,  an  associated 
Y (T^J  value  and  AT^  are  output.  Special  note  is  made  of  any  value  of 
outside  a specified  interval  about  the  assumed  rotation  rate  and  further 
investigations  of  these  areas  are  undertaken  as  necessary.  A mean  AT  value 
is  computed  for  the  well-behaved  portion  of  the  flight  and  output  at  the  end 
of  the  data  pass.  This  AT  value  for  each  lateral  magnetometer  normally  consists 
of  spin  and  processional  motion  and,  therefore,  is  referred  to  a-;  the  rotation 
rate. 

Although  no  specific  mention  as  yet  has  been  made  of  the  axial  magnetometer 
output,  tne  routines  apply  equally  as  well  in  the  case  of  this  magnetometer. 
Computation  of  maximum  and  minimum  values  provide  us  with  a good  approximation 
to  the  cone  angle  of  the  vehicle  axis.  Should  the  output  exhibit  only  the 
slightest  variance  in  successive  maximum  and  minimum  values  respectively,  then 
the  vehic.e  axis  motion  is  due  almost  entirely  to  processional  motion  and  the 
cone  angle  will  remain  nearly  constant. 

If  the  vehicle  axis  study  option  is  desired  in  addition  to  the  lateral 
magnetometers,  supplemental  outputs  include  the  cone  angles  and  precessional 
period  of  the  vehicle  axis.  Examples  of  normal  outputs  are  displayed  in  Figures 

23  A 6 B. 
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Some  of  the  more  common  types  of  problems  with  gyro  systems  include 
situ.itiois  in  winch  i lie  gyro's  roll  output  remain  constant  for  a portion  ol  a spin 
or  lie  y/ro  men su remenl  ■ go  outside  il  calibrated  limits  eithei  with  or  without 
loss  ol  initial  reference.  Lxamples  of  sucli  problems  and  correction  teclmiques 
to i 1 ow  . 


The  vehicles  At) 9. 303-4  and  A10. 312-3  had  onboard  MIDAS  gyro  systems  whose 
gyro  roll  output  remained  constant  for  a portion  of  a spin  (Fig.  24).  To  correct  t: 
problem,  a roll  rate  was  established  for  the  ill  defined  flight  areas.  These 
roll  rates  were  calculated  from  the  roll  data  on  each  spin  that  was  not  influenced 
by  this  malfunction,  and  they  were  substantiated  by  rotation  rates  determined 
from  the  measurements  of  a magnetometer  which  was  mounted  normal  to  the  rocket 
axis.  Utilizing  the  predicted  vehicle  roll  rate,  a lock  was  made  on  the  well 
defined  area  and  continuous  roll  output  resulted. 

The  ACS  system  for  the  vehicle  IC630.02-1A  went  outside  its  calibrated 
limits  (Fig.  25)  on  the  pitch  axis  but  never  lost  its  initial  reference  from  the 
point  of  ancaging.  Converted  reliable  true  pitch  (P  ) data  as  a function  of  time 
was  available  for  various  P.„  values  in  which 

0°  < P < 110° 

25(1°  P < 358°  . 

These  reliable  data  points  were  used  as  samples  for  least  squares  polynomial 
fits  (see  Section  2.3)  in  each  of  the  various  time  intervals  during  which  satura- 
tion occurred.  The  criterion  for  an  acceptable  polynomial  approximation  was  not 
only  a minimum  and  acceptable  RMS  value  between  predicted  and  non-saturated  measure 
P values,  but  also  continuity  between  the  predicted  P.j.  data  in  the  saturation 

areas  and  measured  P...  data  in  the  non-saturated  areas. 

1 

Recovery  procedures  for  vehicle  IC511.21-1A  were  much  more  involved  than 
those  discussed  above.  The  oscillogram  inspection  indicated  that  the  gyro  mal- 
functioned from  250-305  seconds.  To  predict  yaw  data  in  this  interval,  an  in- 
depth  analysis  and  examination  of  the  plots  of  gyro  measurements  suggested 
a modified  Fourier  Series  of  the  form 


A + Bt  + C sincjt  + D cosut 


be  used.  Lxammation  of  the  raw  data  showied  an  acceptable  precessional  period 
to  be  11.2984  seconds  and  accordingly  we  set  the  angular  velocity 

27 r 

M = 1177954  • 

Results  of  computer  runs  produced  a yaw  RMS  of  0.0°.  This  was  deemed  acceptable 
and  the  equation  to  fit  yaw  data  in  degrees  as  a function  of  time  t was  defined  as 

Yaw  (T)  = A + BT  + C sinwT  + D costoT 


where 


T = t-100  and  250  <_  t 308 

A = -14.234  B = -.000  C = -6.354  L)  = -15.838 
and  w is  as  described  above. 


f 

To  predict  pitch  data  in  the  above  mentioned  interval,  plots  suggested  a 
polynomial  approximation  be  used  due  to  the  relative  steepness  of  slope  in 
available  pitch  output.  A thorough  examination  of  plots  indicated  that  a 
structure  in  the  pitch  data  from  250-308  seconds  was  similar  to  measurements  in 
the  portion  from  130-198  seconds  after  launch.  A cubic  equation  yielding  an  RMS 
value  of  .5159°  produced  the  best  polynomial  approximation  to  the  pitch  data 
in  this  well-behaved  area. 

By  the  Weierstrass  approximation  theorem  (Reference  [2],  page  102)  , we 
introduced  the  function 

Ay  (t-271.9)  to [250 , 308 J 

4-0 

and  predic  :ed  pitch  data  in  the  closed  interval  250  < t < 308  by  the  equation 
Pitch  (tj  = AT5  + BT2  + CT  + I)  Ay  , T = t-213  . 
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The  constants  in  the  Ay  function  wore  determined  by  the  available  pitch  data. 
The  values  of  coefficients  determined  bv  a linear  least  squares  estimation 
(Reference  [ 1 j)  were  found  to  be 

A = - . 00034  33064  B = .06455981  C = -3.135508  I)  = -1.314757  . 

Substituting  the  values  of  the  coefficients  and  replacing  T by  t-213,  the 
equation  to  predict  angular  pitch  data  becomes 

1 itcn(t ) = -.  0003433064 ^+.28393068264.“  -77. 32826863tt  -6903. 397458 
with  an  RMS  value  of  .5980°. 

from  the  Roll  Rate  program,  the  vehicle's  roll  rate  was  established  for 
this  gapped  area  and  through  continuity  checks  was  locked  into  the  roll  data. 
All  gyro  corrections  were  incorporated  into  the  aspect  program  to  provide 
continuous  preliminary  attitude  for  the  entire  flight. 

furthermore , a phase  shift  study  between  the  actual  magnetometer  measure- 
ments and  those  simulated  using  attitude  information  was  initiated.  The 
results  of  this  study  showed  an  average  time  shift  between  maximum  predicted 
and  maximum  measured  magnetometer  attack  angles  to  be  .0554  seconds  with  an 
RMS  of  .0048  seconds  for  the  interval  from  launch  to  54  seconds  in  flight.  A 
similar  time  shift  study  from  107  to  272  seconds  showed  the  average  time  shi  ,'t 
between  maximum  attack  angles  to  be  .1694  seconds  with  an  RMS  of  .0182  seconds. 
Phase  shift  displacement  angles  were  also  generated  throughout  the  vehicle 
lifetime,  but  they  did  not  show  a constant  shift.  The  spin  rate  change  varied 
from  .6  rps  to  1.1  rps  from  launch  to  107  seconds  and  then  stabilized  at  1.08 
rps . 


In  an  attempt  to  use  the  times  of  the  maxi mums  and  mini  mums  of  the 
actual  magnetometer  measurements  to  further  refine  the  roll  data,  an  analysis 
was  developed  using  the  yaw  and  pitch  data  at  the  times  of  the  above  mentioned 
magnetometer  measurements,  the  orientation  of  the  magnetometer  and  the  gyro 
reference  notch.  To  implement  this  analysis,  the  yaw  and  pitch  data  was 
merged  wiih  the  times  of  maximums  and  mini  mums  of  the  actual  magnetometer 
data  as  supplied  by  the  detailed  phase  study. 


A software  routine  was  then  prepared  to  calculate  the  discrete  gyro  roll 
using  the  merged  gyro  information  and  theoretical  model  of  the  harth's  magnetic 
field.  The  results  of  tins  program  were  studied  and  continuous  roll  was 
generated  for  this  flight  using  these  calculated  discrete  roll  values. 

In  the  judgement  of  the  Initiator  of  this  problem,  it  was  decided  that 
the  roll  data  calculated  using  the  refinement  technique  of  incorporating  phase 
shift  outputs  compared  most  favorably  with  his  attitude  sensitive  experimental 
results. 

Figures  (26-28)display  the  final  angles  of  attack  of  these  vehicles  using 
the  various  mentioned  correction  procedures. 
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Figure  11 


S33303Q 


28 


50.0  100-0  150.0  ZOO-O  250-0  300.0  350-0  400 

SECONDS  AF  TER  LAUNCH 


31 


Jblib  JO  3 ~1 0 N b 


Figure  18 


CONDS  Rr  It  R 


IHiSF  Si  IFf  IE  f Ufttt  Tif  i!*Slll»eo  #*40  “CNC ttf CO  MSCMtTSHtTt^- OUTPUTS' 


li^uro  23A 


r m 

-i  X 

0-  to 


O to- 

a — » 

*-  t.i 
!3  O 


rp  Os  *-• 

> 0 .»  0 
Jt  -y  .♦ 
OOP 


4-  S v’ 
*»  0 k 
/I  -T  *to 
O Cl  O 


-■  « M »l  P I* 

O "P  k.  Os  tp 

+ J?  & JT  it 

POP  O O C1 


X X 
C.  to 


p 

«r  ' 


0 os  »to 

j »P  r 


O Pj  ru  e»  .♦  > 
rj  o K n 


.£*  ik  O'  P-  t k fO 

*s  0 p.  0 *\j  0 w'  p .r 

-to  -to  -to  J J if  J 

o o p o o p o o o 

* r * ; * r * r 


0 jr  .to  tv.  Is  it  .<  D ki 


INJ  .t  ti 

a>  a»  0 

Jt  J JT 
OPC 


! ! ' I I I I ! 

p (n  k •.  h ^ c~>  _t  nf'jrv*,»»tvpir'orifTi^t,o 

0000k.4f000»i3«c>|f'Cjfsk.  -to  0 fs  0 .j  (Vi  0 0 .to 

J .1  j .*  jr  > j j j jr  ^ .t  .t  .»  j j ^ j ^ ^ j j 

o o n c o n n P o o u P o o u p o ft)  o 6 o cuj  o ci 


jo  os  o. 
O'  44  -t 


«r»  *i  o.  ! 
o j rj  i 


I H U»  Cj  W 


K.  M O K1  > n o ^ <f) 
JOO^tflPPK  O f\J  O I 


Pi  *0  O 

0-  ►*  k. 

! ! 


•i  ^ O *i  N t> 
k ^ s.  k ^ 


ec^r>ircT'iAPo^cvj 
tD  vT)  J)  il'  O \S  w)  «P 


n.M  ^ 
0 s>  0 


-to  4-4  Pj  ( 

0 0 0 i 


1 r C I 
'IAiCi 


O IO 
O O 


O'  m 
t'  >0 


(\J-iKBr,Jr‘N.vP^fr>ff' 
vD  vD  IM^  vD  IT  t''  If'  O l/-. 


Os 


10 

*- 

«? 

or 


C- 

**> 


k . ki  r\j 
iD  Jf  K 


^ *»  i 

• w • 
-*  -*  *» 


O'  r»  ^ 
O'  0 too 
it  «c 


IMT  o o ^ 

-i  © ^ ^ IT 

I 00  *4  O Pi  CP  O 
I • *•  • >•  • N 

J1  0 ^ ^ f 


0 O'  f\i  IP  m if-v  ift  s-<  .o 
w v k.  o r ■ vfi  44  in 
r.»ir-HO\^J  < T 
O'O'CPO'rT'tptnOsaO 


-to  <r>  0 
PJ  -4  ** 

eo  #s  p. 


*r>  Os  44  Jr  p. 


o CT  (\j  i 

K.  «5  O'  f*i  P P 
Os  ov  0 »v  k.  p- 


0 c . 
o o»  < 
O vP  I 


. 0 »s 
10  10 


ww«mwmMif)wnnmK)«r> 


o-  #»  ff1  ih  P1  frJ  o # 
**  rv  n*  0 -to  0 ■*-« 

♦-  l\J  fU  M IP  M -to 
fj  M CVJ  f\J  »J  M IV- 

* t*  * r 


a 

=>  lu 

CO  *J 


IT  i 


« $ 


0 >0  0 
cr  <y*  tr* 


i K.  Os  CtJ 

I r>  CV  "> 

I w ^ •o 
k • «• 
0 0 0 
it  cr  r 


« 

r 

I i 


00*0 
k.  (T*  rj 

k tr  u 


o o o 

(D  0 *t 
* %£>  & 


K.  iw  f\j  iMP  J) 
»<  C k I.  K ifl 


rv  oj  pj  Oo  rv  ou 

• {•  •(••■• 


!.*«  W O J lj 

r*  .to  **  +*  ® c> 

in  □ (j  o d 

• • • »•  • k 

OS.  \Ck\fiK. 

^ <y  <y  j'  cs  y 

! i i 


»<  O k O 
oo  je«s 
-t  0 o ® ^ 


I j » I 

.T  *f>  « (Vl  «)  ir  m n»  K 
a w o ^ p ^ O'  ov  0 
iT'U'lftU'UMMniMT 


I 


i in  O o n i it 

• 4-4  »0  Pj  Pi  Ok  O' 
I n © V P ^ w 
• • • » • » • 

K \C  vD  S.  O K 
cr  O'  O'  O'  CT  O'  O' 


10  cr  CVJ 

— «T.  > 

*d  0 »n 

f J fVI  0W 


• r 


I 

rvi  -*■  «n 
m k.  o 

0 10  K. 
fVJ  f*J  PJ 


n n n r>  n 


rr>  n 


i 

.♦  w» 
O fs 
tO 

0 0 

• * 
*n  K> 


.♦  cvi  *4  ip  •*«  ci  ru  —4  p.  . 
k «*<  O O.'  .0  ♦<»»-•.»  K.  , 

000  1^  00  00010' 


• V • • 


(t  «t>  m •<  e iVi  *4 

fK  fT  k •<  O'  K>  n 

if  iP  k k 0 »w  «p 

(\j  lij  M tv  ry  (\i  (vj 

• t 


; 

k-  a pj  f.  0»4fi.  PS  ^ O'  o M 1 

LO  0 iti  n © r j m pj  ns  ^ if, 

k-  n»  ►»  k.  «o  n»  k tv  an  k.  <n  k.  i 

OJMtVfvrJfVtVtNjfVMrVfyfVil 


k-  0 

o'  cy* 

I 


I 


cr  io 

Pi  ec 


•o  O'  -r  !»■»  cvj  O'  0 
OiCtU® 
•4  W>  0 O' 

• •»  • W 


H t *C 


(\j  in  tv 
® ao  o 
Wk  iD 

k-  0 vS 
O'  cr  O' 


r>-  0 *4 
^ ^ «> 

O'  (NJ  ^ 


I 


O'  K)  . 

tv«  o i 
0 0 0 1 
• 0 • 
k.  k.  0 
O'  O'  O' 


O'  O' 


T 


m 10 

H 0 
• » 

k.  0 
O'  O' 


iD  *4  Ti  ♦< 

0 **«  0 0 l 

k ct>  cw  0 . 


O'  O'  o ( 

Pj  f\J  ■*>  I 


I -4  *-  «4  PJ  f>J  l 

' K*  IO  »*>  l*>  ■ 


I ^ W P J 4 J . 

i w io  w t*>  r>  Hi  m I 


>0001 
l W W f*)  I 


• h 

• k.  i 


00  O' 

0 »P 
*P  0 


4rk#iTtiiiMy'i\iin0^Dct 
♦4  10  o k « ►»  w o ^ o ru  0 ir 

•44fOO'«500-fV0^H0«0»S.l 


rw  ^ 
«t>  0 


k.  Ok 

O'  C0 


ru  fo 
0 0 
^4  J- 


O'«)0rsfv0^40^;0fkor^ 

♦••••©•©•OoOofc 

0k,o*000e0tc»rkK.0OkOv. 
O'  O'  O'  O'  cpO'O'O'O'CT'C’O'C' 


IP  OP  I 

I 


li)  *4  fJ  fo  u 

0 J 

0 •< 

4 fli  »4  ^ (V 

m^»in  i 

•4  0 

0 . 

(V  O ffl  X 

lokin  « 

•-*  P 

0 r 

pj  in  k e 

#)  O'  cr  O'  n 

o c 

> oo  « ^ ^ ru 

pp«ip  f 

> -y  J 

•O  J 

^ o 

op  ip  cr* 
^00 
Ok  0 *4 
> > ^ 


•4  »4  0 
»<  k / 
O 0 0 


k.  ♦*  0 0 O'  .y 

•P  IT-  -4  O C“>  **-> 
O'  IP  «4  ♦-  er»  C5 
Pj  «4  ^4  <VJ  O'  W 


»-4 


10  O'  10 
O'.  KS  0 
1 Pj  Pi  Pj 
I P>  Pi  PJ 


10  PJ  k. 
PJ  tC  O 

io  pvi  -r 
Pi  Pi  Pj 


'*'  X k u'  k k 


P ^ ^ P If  ^ 
Pi  Pi  Pi  At  Pt  m 


o-  U* 

T>  Q 
O 

o 

UJ 

0- 

-T  U 


«c  o 

O' 

Pj  4-1 

iii 

CM 

0 

*■4 

K 0 

m 

4-4  ^ Is 

0 

0k 

^4 

O' 

0 

PJ 

O 0 

ip 

0 0 4-4 

> 

J' 

0 

1^  J* 

op 

lO  0 

Pi  Ik 

o 

O'  O' 

Vp 

IP 

IP 

PJ 

-P  H 

O' 

0 JO  IP 

0s 

PJ 

10 

*4 

o 

r> 

0 tr 

i 

s 

k.  0 k. 

■t 

CS 

o 

—4  CO 

0 

«4  44 

Pi 

° 't 

CU 

J- 

0s 

IP 

^ Pi 

0 

d k O 

«-t 

IP 

Ov 

Pi 

«r 

Pj 

0 O 

M 

0 k Pj 

0- 

cr 

Pj 

0 O' 

0 

0 Os 

O O' 

• 0 

-* 

• 

o*  e 

» • 

T 

0 

f£ 

• 

0 

«f)  k» 
• • 

Ik 

to 

Ik  Os  0s 

• o • 

tk 

0 

k. 

0 

► 

0 

0 

to 

k. 

• 

0 0 

0 0 0 

• I • 

IS 

0 

0 

» 

0 0 

• 0 

0 

0 0 

^ IP 

«* 

ip  ip 

IP 

| 

op 

r> 

IP 

IP  IP 

I 

OP 

Ip  op  IP 

ir> 

op 

IP 

IP 

IP 

ip 

K)  |p 

ip 

IP  IP  IP 

1 

IP 

OP 

f 

3 

- -4 

IP 

K»  IP 

PJ  Ip 

Ov 

.t  «r 

Os 

Pi 

i 

r» 

0 

4*4  k- 

0»n 

10  IP  o 

r> 

O' 

04 

k* 

cr 

0 0 

0 

k.  o Os 

ft 

ip 

J, 

1 

-4  O' 

4? 

cr  0 

m 

0 O' 

-r 

r-t 

0 

O' 

cj  rr> 

IP 

Ip  of'.  Cj 

cr 

Ik 

cv 

0 

4*4 

tr 

ro  is 

0 

Pj  a*  0 

O 

sr 

.p 

Ok 

ft 

CO  o 

!y  00 

-r 

0 l*' 

10 

0 

0 

0 

0 0 0 

0 0 0s 

0 

0 

tv 

k. 

0s 

0 

Is  0s 

Is 

kkk 

Os 

Os 

Ov 

Is  k. 

•P 

0 an 

pj  pj 

Pi 

rvj  pj 

pi 

PJ 

Oil 

Pi 

Pi  Pi  PI 

Pj  Pj  Pj 

pj 

Pi 

Pi 

PJ 

fll 

Pi 

Pj  pj 

r 

Vi 

pj  rv  pj 

cv 

Pj 

ft 

Pi  rv 

Pj 

PJ  Pj 

Oil  P.' 
rt  tr 
K k. 


I 


I 


P p ^ k (V  k k t «■  «4  O'  fi  0 0 0 10  0»h  cD  (\j  ifc  0 k-  (\j  «<-  c*  r»  J O0« 
^4  pj  t»  0 a/»4P  fjikO4f,j0Kk  T00»40«-4CT'C»0O'0O'o-»«t;} 


O O'  i0  O'  O'  -t  r*>  0 

•4  in  i p .t  O'  pj  i*5  O' 


O'  O'  »4  ^ .t  0 m X in 
.t  tr  k 4 if  o m o -j 


or 

*J 

P 

<r 

cr 

O' 

•4 

0 

0 

o* 

Ip 

IP 

0 

tr 

0v 

0 

Ov 

K>  0 

IP  0 

OP 

0 

«?  0 

0 

J 

0 

00 

0 

0 

Os  0$ 

4» 

ru 

0 

P 

-to  to  pj 

00 

Ov 

CJ 

0 

O' 

0 O'  Os 

tv 

l*. 

IP 

w 

• 

• 

• 

s 

• 

s 

• 

■ 

• 

• 

■ 

• 

• 

• 

» 

• 

• 

• 

• 

• 

• 

0 

• 

• 

• 

• 

• 

to 

• 

• 

• 

0 

• 

■ 

• i 

• 

• 

• 

0 

• 

I 

• 

0 • 

• 

r 

r 

O 

n 

o 

r-» 

n 

*4 

f) 

*4 

o 

IP 

44 

CP 

♦4 

rp 

CP 

4-4 

IP 

IP 

•44 

•4 

r> 

H 

4-4 

CP 

o 

4-4 

-4 

rp 

rP 

4-4 

O* 

rp 

44 

44 

44 

rp 

4*4 

*4  «J 

1 4-4 

o 

*4 

44 

Pj 

44 

r> 

O CP 

44 

UJ 

cr> 

r— 

rrt 

r 

cr* 

o 

IP 

IP 

o 

IP 

o 

«T 

o 

ip 

o 

C> 

O *"• 

O' 

n 

ip 

CP 

CP 

CP 

o 

n 

CP 

CP 

o 

CP 

cr* 

rp 

CP 

e> 

CP 

CP 

o OP  o CP 

o 

O 

o 

CP 

rr 

CP 

Ct>  rp 

cto 

o 

«4 

»4 

4*4 

^4 

•* 

*4 

*" 

•* 

** 

*4 

*** 

*-4 

*"* 

••4 

4-4 

44 

4*4 

4-4 

•f 

44 

44 

4-4 

44 

•4 

4*4 

*■* 

•" 

to4 

44 

•4 

•4 

*4  4H 

1 44 

44 

*4 

4-4 

%* 

44 

44  H 

1 

7 

k 

♦ 

vs 

0. 

•4 

0 

Os 

IP 

k 

O' 

Pi 

k > 

1 

*s 

ip 

p. 

P 

.K 

k. 

A' 

CP 

0 

rr- 

Os 

0s 

.rt 

- 

[ 

0 

Pi 

0 

if* 

tf  |fl 

1 44 

j 

-4 

t » 

ro 

Os 

an 

>* 

.♦ 

T 

O' 

w- 

l s 

ar 

v. 

4» 



. * 

f 

IP 

*r 

10 

r 

s. 

cr 

to 

ip 

IP 

_ 

o* 

.t 

o 

Os 

mt 

p 

■ m 

OA 

r-> 

* 

♦ 

P 

O' 

0 

IP 

rr 

»P  O- 

► 

O' 

0 

to 

. 

ro 

k 

• » tr 

ir 

k» 

CT 

■9 

0 

«C 

«•• 

CP 

tr 

i/ 

c< 

CP 

IP 

0 

cr. 

44 

»*> 

if 

O' 

•* 

0 

0 O' 

Pi 

IP 

Os 

o 

OP 

IP 

0 

**• 

-to 

tr 

O'  PI  0 

P. 

c * 

tr 

0 

«o 

4*4 

4 is 

r* 

PJ 

rr 

r 

O* 

r 

f. 

r* 

IP 

— 

*1 

— 

-1 

P/ 

»*  • 

p • 

P* 

•p 

tr 

tP 

•p 

♦ 

to 

to 

ip 

IP 

ip 

IP 

.P 

>n 

»r 

Ov 

Is 

Os 

0s 

«r 

.f 

0 

r P m 

t2 

CP 

rs 

ft 

r» 

si 

2TH 

tv 

PI 

IP  OP  OP 

I 


ppppppppp^ppMpPKipK>PPPPP»«P«MPPfOPl^PPj  .t  jf  > .t  .♦ 


40 


GEHER&TEO  OUTPUT  HEISUREO  OUTPUT  ®xiSE  TIME 

-M»*-  rilE  — **4  Jt  E OEtnr  — 3PIN  R»TE HRR-TINE *NltE 3*trf«T — SPIK  *«Te— -S-TEf SUFI 


Figure  23B 


© 


0 


i k 

qr\  ►«.  vP  k. 
■k  -*  4r  * 
(j  cm  o 


♦.i  r»  X*  -* 

f'  w • vT 


**■'  «x>  to--  ^ ir>  *4  m N.  # 

O N-  ^ Q)  4 ifl  M 

c>  e>  o o o o o o < 


-*  - r - if  -t  J J J 

o ci  u o o o u n p 


I 

<Yj  .4  *.  rvj  pf»  .♦  .y  k.  tr> 
tr>  O'  rv-  «r>  k O'  «•  O'  tf' 

t ^ ^ ;r 

OOf}D(tiUOO<> 


^^1  oi  1 ’ 


4«  «4  «4  k*  «4 


O'  *-  ♦*  ri  l/>  ^ 

» k »■  p-  i 


+ ,n  rsi  ,-r  jr  **- 


> in  iH  . 
k rl  i 

1 V ' 


**  K pn  ^to  »*•.  «• 
> >•  • 

H>  »>- 

i ♦«  • « 


. 

■w  N A'  ifi  <h 
‘ OJ  O ♦<  ID  W iJD  I 
I 04  ru  Oi  Vi  1 

> « O O O0  O i 

1 *•  • k • *•  • 

I ^4  «4  Y *4  +*  *4  ' 


MN  Ki  0,m  (i\  O ft!  ,n  ^ W ^ > M 


* o © i 

k 


*Ur>di<T'<f*»4»v,^«3 
PW®IS.»**fW®P'.«K- 
• O'  0 ^ 0 fP  ^ O'  O*  ^ 


CM  5J 
m <t> 
O*  O' 


e»  > < 

«n  » i 

O*  O'  < 


I « o 
i rvj  pt> 
«r>  «o 
, O'  O' 

't* 


I 

,#-  *^  cu  -j  di 

1 

I lO  U>  tp  IMM^ 

!(••*•  • » • 
I (un 

i © © © ® <5  © 
I44M«4«4^4H 


O'  ® 
iD  IT 
-f 

• » 

i «4(M  I 

■ O'  O'  < 


■*  K 


© c*  O'  tr  -k  p'.  O'  © © 

K K K K h K I 


^OOO0O0OC»I 


*toW»4«|»«|totop*f»»»#^' 

-'H«4*4«4*444^4*4k4to-«toM«>4«4«4©W44«4«4 


:t":t 

I 4t«4  4c 


) 1 O K I 


t0^ 


PO  4>  in  *4  i . 

KiN0ibfr0rurf>H,  . 

«0©«D©©©0'®0'< 

O'O'O'O'O'O'O'O'O'C 


I 

00^K000nH-a«r4«tr0VN0U'0oN.ocriNjit0. 

0®»pKWl\i*c*C4'OlV*C>DCr'  ff>lMrW4Mr^*CfVOfW®l 

0'0'©0*4<Ni<ur>PP>.?.yininininU>P'*Pto.«>©0'0'©«*«f««4 

«Vi«M»f»*opr»>oo>wo>R*>»r)#o»r>rofpfow>0>»0iorT>»0-t^jf^ 

«4«*4*v«44«44««44<<PC<*(«<^*«4CH«HH«CvC44«4«H*4«k«« 

II 

-k'£l«bO'»®m<>-»-»fto.O'P'».OPV*£ru»4©r\J«»«in 

^N®e«^(uo0tOir®<4O®ir)tf'^nj«4®9'4k^tft 


i ^ o jr  »<  ftj  •? 

iQfyqjo^e 

I N N W N (\j  ^ 

loeedoo 


•*  v£) 

O'  Pv 


-y  ru 

O'  ® 


• ->  in 
> O'  o* 

■r- 


K »♦)  W 1 


O(V^n>«40M«M(M 
^ru«4OO'«0K.4>tn. 
nwwff)cjrij(\jrv(vji 

'0*0000000000 
i«^«4<yr\jc\JCM<v<Mfv<\J<)J 


°1 


»D^it)SKOMrt 

■O  y m cv  -4  0 o O' 


|,r«,-l,Hv4,W4.»4«4^.Mv4««4'«4«4«-««4*40O 

iO0o600o0OoO0O0ei 


irs  > 
p*.  tv 
5 


O'  »4  vb  *0 
Pw  c-.  *.  O' 
N-  B ><•  K 
O'  O'  O’  O' 
»•  • f • 


e (Mp  w m 
© © © © O •* 
© K.  © © Pw  © 
O'  ^ O'  (T  ^ jj* 

k • p • • • 


i. 


> © «> 

' O'  O' 


/“»  '•*4  1 
-4  TM  l 

• «*  < 

O O'  ( 


O'  O' 

1 r • 


O'  l/' 
it>  fO  ">  «> 
«44  «D  l/'  O' 
» • • • 

v£>  K.  Pw 

P0  f0  Pf)  »0 


*40  ® 

fvj  y y o f\c 
(T  4 ® C PO  N 

• O'  O'  e>  «m  *4  f 

® ww  yy  ^ • 


* .y  d>  . 

I vP  O l 
PO  \P  . 


P0P0P0P0P0P ^ ^ ^ -t  ^ ^ 


-4  ^ 
p0  in 
4*4  P*. 

i/  in 
-y  -y 


ru  «4 
o»  rvj 
ro  4-4 

to  • 

vD  p^ 

«t  y 


l 

P0  O' 
-r  in 
«n  «o 

O'  O' 


•»•••*•* 

cR  0 vi)  IP  y (P  i 
if*  1*1  ^ «r  O O'  t 


®®ir«r0B(to 

O'  O'  O'  O'  O'  OP  O'  O' 

r • > • r • r 


O'  O' 
O'  <T 


O'  vD  *0  ^ 
P*  fM  Jt  « 
«P  Oj  CT  sD 
► • 9 • 

K «?  «5  O' 

.y  jy  ^ j- 


n.4-4vf>in«©*£in**»vt> 
o'4fM®oyyvCoj 
PJON  PHVOifl  JN 
Pto^toi*  • • to 

rt»44^nj**i^y^i^vfc 

u>ininmu>inu>in 

to)  v4  4*  *4  Jl  *4 


4 * 0 4* 

K K0  B 


to.  .to  4)  i£  IP  to. 

I O'  Pto.  in  »n  W O I 
IP  IMP  IMP  4 

*•»•»• 
toiw  y in  4>  to. 

O oo  «t>  w o c> 


n ••*. * 

•n  ♦■  ( 

-y  -r  ■ 


i y ® i 

• — . 00  I 

y PR  i 


1 k ru  ® Pto 

i N n n 

PO  no  po  nj 


in 
v I 

Jr 


^ tor«  . 4 0<  *0  .♦  ito  (p  k.  ® fto 

<C®(J’ff'ft(P(T'0'J’(P(*ff* 

4»4  «4  «4 


© ® © O © 0 © 


m to*  nj  to)  .pi  «r 
• l#.  .♦  PO  fM  m*  CT 
M»ufvfu^(\jNnj*H 
Vto^toPtototoL 
n w <y  m y os  to  k.  4i 

0oo  oo  o 0a  k 

PPMrylUNMftPMfy 


O © 0 © 


» 

1 

I 

I 


9 


i 

I 


I 

i 


41 


-tut  4W?RJ  JE-TIHc  SHIFT  IS w5VP5  -MITM-Vt-RHS-flF rflOEV 


HNGl.f.  or  RT  TRCK 


L HUNCH 


I 


RhFP.RhNCFS 

1.  Daniel  C.,  and  Wood,  F\,  Fitting  liquations  to  Data,  Wiley  and  Sons 

Inc.,  1971. 

2.  Macon,  N . , Numerical  Analysis,  Wiley  and  Sons  Inc.,  1963. 

5.  Roxborough  et  al , Procedures  for  the  Determination  of  the  Attitude  of 
a Rocket  from  Gyroscopic  Data,  Boston  College,  Final  Report 
AFCRL-72-0740 , prepared  under  Contract  Number  F19628-70-C-0017, 
September  1972. 


A 7 


